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This  paper  presents  an  electrochemical  model  for  simulation  and  evaluation  of  the  performance  of  proton 
exchange  membrane  (PEM)  fuel  cell.  The  results  of  the  model  are  used  to  predict  the  efficiency  and  power 
of  the  fuel  cell  as  a  function  of  operational  parameters  of  the  cell,  like  temperature,  partial  pressures  and 
membrane  humidity.  The  influence  of  temperature  on  fuel  cell’s  characteristics  is  more  pronounced 
than  the  influence  of  partial  pressures  and  membrane  humidity.  The  effect  of  platinum  loading  on  cell 
performance  is  examined  with  Pt  loadings  of  0.1 8, 0.38  and  0.4  mg  cm-2.  The  kinetic  parameters  (electron 
transfer  coefficient,  exchange  current  density)  are  found  to  be  platinum  loading  dependent. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Within  such  perspective  PEM  fuel  cell  systems  have  been  show¬ 
ing  up  as  a  promising  alternative  due  their  high  efficiency  and  low 
impact  to  the  environment.  The  theoretical  approach  with  mathe¬ 
matical  simulation  is  very  effective  for  designing  and  analyzing  the 
performance  of  the  PEM  fuel  cell. 

Some  work  has  been  reported  in  the  literature  on  steady-state 
fuel-cell  modeling  (e.g.  [1-4],  as  well  as  dynamic  modeling  [5-7]). 
These  studies  are  mostly  based  on  the  empirical  equations  and/or 
the  electrochemical  reactions  inside  the  fuel  cell.  Also  we  can  find 
several  researches  regarding  influences  of  certain  operating  param¬ 
eters  on  fuel  cell  characteristics  [8].  There  are  many  other  models 
have  been  developed  and  reported  in  literature,  but  they  did  not 
focus  on  the  impact  of  individual  operating  parameters  on  fuel  cell 
output  characteristics. 

This  paper  presents  a  simple  model  of  a  PEM  fuel  cell  that  can 
be  used  to  analyze  the  impact  of  individual  fuel  cell’s  operating 
parameters  on  cell’s  performance.  Using  the  present  model,  it  is 
also  possible  to  determine  the  activation  loss  parameters  of  the  PEM 
fuel  cell.  The  model  is  well  adapted  for  PEM  cell  and  it  incorporates 
the  essential  physical  and  electrochemical  processes  that  happen 
in  cell  along  its  operation. 
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2.  Experimental 

2.  I .  Electrode  preparation 

For  the  electrode  preparation,  a  substrate  made  from  car¬ 
bon  cloth  “A”  with  thickness  350  |jim,  platinum-supported  carbon 
10  wt%  Pt/C,  activated  carbon  black,  and  Teflon  (PTFE)  60  wt%  have 
been  used.  Two  different  platinum  loadings  have  been  used  in 
preparation  of  the  electrodes.  One  electrode  was  with  platinum 
loading  0.1 8  mg  Pt  cm-2 ;  while  the  second  was  with  platinum  load¬ 
ing  0.38  mgPtcm-2.  The  composition  of  catalyst  layer  was  70  wt% 
Pt/C  and  30  wt%  Teflon  (PTFE).  Detailed  electrode  preparation  can 
be  found  in  Ref.  [9]. 

2.2.  Operation  of  single  cell  in  fuel  cell  test  apparatus 

During  operation  of  PEM  fuel  cells,  the  following  processes  take 
place  within  the  electrode:  (i)  the  reactant  gases  diffuse  through  the 
porous  backing  layer;  (ii)  at  the  gas-electrolyte  interface,  the  gases 
dissolved  and  then  diffuse  to  the  electrolyte-electrode  interface; 
(iii)  electrocatalytic  reaction  inside  the  catalyst  layer  precedes  the 
gas  adsorption  at  the  electrode  surface;  (iv)  ionic  transport  occurs  in 
the  electrolyte,  but  electronic  transport  takes  place  in  the  electrode. 

The  oxygen  and  hydrogen  were  passed  through  humidifiers 
before  being  fed  into  the  cell  cells.  The  hydrogen  fed  into  the  anode 
at  a  flow  rate  of  140  ml  min-1  and  1  atm.  The  oxygen  entered  the 
fuel  cell  through  the  cathode  at  a  flow  rate  of  380  ml  min-1  and 
2  atm.  The  electrons  generated  from  the  anode  were  connected  to 
a  digital  multimeter,  with  an  external  variable  resistance  to  mea¬ 
sure  the  current  and  voltage  produced  by  the  cell.  The  electric 
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Nomenclature 

V activation  activation  voltage  drop,  V 

Vohmic  ohmic  voltage  drop,  V 

V concentration  concentration  voltage  drop,  V 

E  fuel  cell’s  reversible  voltage,  V 

F  Faraday  constant,  C 

i  current  density,  A  cm-2 

io  exchange  current  density,  A  cm-2 

z0c  cathode  exchange  current  density,  A  cm-2 

kc  factor  related  to  reaction  speed,  A  cm-2  C-1 

iL  limiting  current  density,  A  cm-2 

ld  diffusion  layer  thickness,  cm 

lm  membrane  thickness,  cm 

Pcell  fuel  cell’s  power  density,  W  cm-2 

Ph2  partial  pressure  of  hydrogen,  atm 

Pq2  partial  pressure  of  hydrogen,  atm 

R  universal  gas  constant,  8.314 J  K-1  mol-1 

rion  ionic  resistance,  £2  cm2 

rei  electronic  resistance,  £2  cm 2 

T  temperature,  K 

V  voltage,  V 

Greek  letters 

\//  electron  transfer  coefficient 

O  membrane  humidity,  % 

X  membrane  water  content 

P  symmetry  factor 

crd  diffusion  layer  electronic  conductivity,  A  V-1  cm-1 


conductivity  was  measured  by  using  resistivity  meter  (Loresta- 
GP  MCP-T600).  The  specific  resistance  of  the  gas  diffusion  layer 
composed  of  70wt%  Pt/C  and  30wt%  PTFE  was  measured  to  be 
0.21  £2  cm.  The  operating  parameters  of  the  fuel  cells  are  listed  in 
Table  1. 


Table  1 

Operational  parameters  of  the  fuel  cell. 


Parameters 

Value 

Temperature,  T 

298  K 

Electrode  area 

25  cm2 

Membrane  thickness,  lm 

40  fxm 

Diffusion  layer  thickness,  ld 

350  (xm 

Po2 

2  atm 

Ph2 

1  atm 

Diffusion  layer  electronic  conductivity,  crd 

5  cm-1 

3.2.  Model  formulation 

The  output  voltage  of  a  single  cell  can  be  defined  as  the  result  of 
the  following  [10,5]: 


V cell  ~  E  ^ activation  V ohmic  ^ concentration  (4) 

In  the  above  equation,  E  is  the  thermodynamic  potential  of  the 
cell  and  it  represents  its  reversible  voltage;  Vactivation  is  the  voltage 
drop  due  to  the  activation  of  the  anode  and  cathode,  a  measure 
of  the  voltage  drop  associated  with  the  electrodes;  Vohmic  is  the 
ohmic  voltage  drop,  a  measure  of  the  ohmic  voltage  drop  result¬ 
ing  from  the  resistances  of  the  conduction  of  protons  through  the 
solid  electrolyte  and  the  electrons  through  its  path;  and  Vconcentration 
represents  the  voltage  drop  resulting  from  the  reduction  in  concen¬ 
tration  of  the  reactants  gases.  The  fuel  cell’s  output  power  density 
is  given  by: 

Pcell  —  ^ celd  (5) 

Each  one  of  the  terms  of  Eq.  (4)  is  discussed  and  modeled  sepa¬ 
rately  in  the  subsections  that  follow. 

3.2.1.  Cell  reversible  voltage 

The  fuel  cell’s  output  voltage  is  determined  by  cell’s  reversible 
voltage  that  arises  from  potential  difference  produced  by  chemical 
reaction  and  several  voltage  losses  that  occur  inside  a  cell.  The  fuel 
cell’s  reversible  voltage  is  a  function  of  temperature  and  partial 


3.  PEM  fuel  cell  modeling 

3.1.  Basic  fuel  cell  operation 

The  PEM  fuel  cell  consists  of  membrane,  which  is  ionic  conductor 
and  through  which  hydrogen  ions  diffuse  from  anode  to  cathode. 
Beside  membrane,  fuel  cell  has  two  catalyst  layers,  made  of  plat¬ 
inum,  gas  diffusion  layers,  gas  channels  and  current  collectors.  The 
schematic  of  a  PEM  fuel  cell  is  shown  in  Fig.  1. 

The  hydrogen  flows  through  gas  channels  on  anode  side  of  a  fuel 
cell.  It  diffuse  through  diffusion  layers  catalyst  where  it  oxidizes 
according  to  the  following  reaction: 

H2^  2H+  +  2e“  (1) 

The  protons  are  conducted  through  membrane,  which  is  insu¬ 
lator  for  electrons.  The  electrons  pass  through  current  collectors 
and  external  electric  circuit  to  cathode.  On  cathode  side  of  a  fuel 
cell,  protons,  electrons  and  oxygen  that  flows  through  cathode  gas 
channels,  react  according  to  the  following  reaction: 

2H+  +  2e~  +  j02  -*■  H20  +  heat  (2) 

Therefore  overall  chemical  reaction  occurring  inside  PEM  fuel 
cell  is 


Anode  Cathode 


Current  collector 


Gas  channels 


Gas  diffusion  layer 
Catalyst 


Membrane 


H2  +  j02  -»■  H20  +  heat  +  electrical  energy 


(3) 


Fig.  1.  Schematic  of  a  PEM  fuel  cell. 
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pressures  of  reactants  and  product  as  is  shown  in  the  following 
equation: 

a  c  a*s  RT 

E  =  +  2 f  -  W)  +  2F  [ln(pw2 }  +  2  ln(P°2  )1  (6) 

where  AG  is  the  change  in  the  free  Gibbs  energy;  Fis  the  constant  of 
Faraday;  A5  is  the  change  of  the  entropy;  R  is  the  universal  constant 
of  the  gases ;  while  Ph2  and  Po2  are  the  partial  pressures  of  hydrogen 
and  oxygen,  respectively.  The  variable  T  denotes  the  cell  operation 
temperature  and  Tre/the  reference  temperature.  Using  the  standard 
pressure  and  temperature  (SPT)  values  for  AG,  A5  and  Tref,  Eq.  (6) 
can  be  simplified  to  [10]: 

E  =  1.229-0.85  x  10“3(T  -  298.15)  +  4.31  x  10“5T[ln(PH2 ) 

+  pn(P02)]  (7) 


3.2.2.  Activation  voltage  drop 

The  activation  polarization  is  related  to  the  energy  barrier  that 
must  be  overcome  to  initiate  a  chemical  reaction  between  reac¬ 
tants.  At  low  current  draw,  the  electron  transfer  rate  is  slow  and  a 
portion  of  the  electrode  voltage  is  lost  in  order  to  compensate  for 
the  lack  of  electro-catalytic  activity.  The  expression  for  activation 
losses  is  given  by: 


which  is  influenced  by  the  membrane  preparation  [11].  The  param¬ 
eter  A  can  be  related  to  membrane  relative  humidity  0  [12]  by  the 
following  expression: 

A  =  0.043  +  17.810  -  39.8502  +  3603  (12) 

The  electronic  resistance  can  be  written  as: 

rei  —  (13) 

ad 

where  ld  is  diffusion  layer  thickness  and  crd  is  diffusion  layer  elec¬ 
tronic  conductivity. 

3.2.4.  Concentration  or  mass  transport  voltage  drop 

The  mass  transport  or  concentration  polarization  results  when 
the  electrode  reactions  are  hindered  by  mass  transfer  effects.  In  this 
region,  the  reactants  become  consumed  at  greater  rates  than  they 
can  be  supplied  while  the  product  accumulates  at  a  greater  rate 
than  it  can  be  removed.  Ultimately  these  effects  inhibit  further  reac¬ 
tion  altogether  and  the  cell  voltage  drops  to  zero.  The  expression 
for  fuel  cell’s  concentration  losses  is  given  by: 

Wentration  =  — 5ln(l  — ^-)  (14) 

ii  represents  limiting  current  density.  This  parameter  describes 
maximum  current  density  that  can  flow  through  electrode. 


RT 

2xj/F 


(8)  4.  Results  and  discussion 


\l/  is  electron  transfer  coefficient,  and  is  unit  less.  This  value 
describes  the  proportion  of  the  electrical  energy  applied  that  is  har¬ 
nessed  in  changing  the  rate  of  an  electrochemical  reaction.  It  is  this 
value  that  differs  from  one  material  to  another,  i  represents  cell’s 
current  density,  whereas  i0  is  exchange  current  density.  i0  is  the 
value  on  the  Tafel  plot  when  the  current  begins  to  move  away  from 
zero. 

Because  of  higher  anode  exchange  current  density,  cathode  acti¬ 
vation  losses  are  significantly  higher  so  anode  activation  losses 
are  negligible.  The  value  of  cathode  exchange  current  density  also 
depends  on  operating  parameters  what  is  shown  by: 


io  c  —  2  Fkc 


exp 


2.46  /3F\ 
RT  J 


(9) 


P  and  kc  are,  respectively,  symmetry  factor  and  factor  related  to 
reaction  speed. 


The  concentration  polarization  results  from  restrictions  to  the 
transport  of  the  fuel  gases  to  the  reaction  sites.  This  usually  occurs 
at  high  current  because  the  forming  of  product  water  and  excess 
humidification  blocks  the  reaction  sites.  This  polarization  is  also 
affected  by  the  physical  restriction  of  the  transfer  of  oxygen  to  the 
reaction  sites  on  the  cathode  side  of  the  fuel  cell.  The  concentra¬ 
tion  polarization  can  be  reduced  by  using  thinner  electrodes  which 
shortens  the  path  of  the  gas  to  the  sites  [13,14].  Jordan  et  al.  [14] 
have  observed  a  dramatic  change  in  slope  of  the  voltage  versus 
current  density  plot  using  air  oxidant.  Such  a  change,  indicative  of 
a  diffusion-limited  reaction,  was  not  so  apparent  when  pure  oxy¬ 
gen  is  used  as  the  oxidant.  The  same  behavior  was  also  observed 
by  other  researchers  [15].  This  is  consistent  with  the  experimen¬ 
tal  data  showed  in  Fig.  2.  Flence,  in  parameter  estimation  we  have 
neglected  the  effect  of  concentration  polarization  drop. 

4.1.  Determination  of  the  activation  loss  parameters 


3.2.3.  Ohmic  voltage  drop 

The  ohmic  voltage  drop  (or  “Ohmic  polarization”)  occurs  due 
to  resistive  losses  in  the  cell.  These  resistive  losses  occur  within 
the  electrolyte  (ionic),  in  the  electrodes  (electronic  and  ionic),  and 
in  the  terminal  connections  in  the  cell  (electronic).  Since  the  stack 
plates  and  electrolyte  obey  Ohm’s  law,  the  amount  of  voltage  lost 
in  order  to  force  conduction  varies  mostly  linear  throughout  this 
region.  This  is  the  working  region  of  the  fuel  cell. 


^ohmic  ~  Krion  +  rel) 


(10) 


We  can  calculate  ionic  resistance  using  the  following  expression 

[10]: 


181.6 


—  In 


1  +0.03i  +  0.062(3L)2i2.5 


(X  -  0.634  -  3 !)  exp  [4.18  (1^) 


(11) 


The  parameter  lm  is  membrane  thickness.  The  ionic  resistance 
represents  PEM  fuel  cell’s  membrane  resistance  to  transit  of  hydro¬ 
gen  ions.  This  resistance  strongly  depends  on  membrane  water 
content.  The  membrane  water  content  is  described  by  parameter  A 


The  experimental  data  of  the  cell  voltage  versus  current  den¬ 
sity  (Fig.  2)  for  the  two  fabricated  electrodes  and  the  commercial 
electrode  E-TECK  with  Pt  loadings  0.18,  0.38  and  0.4  mg  cm-2, 
respectively,  are  fitted  to  the  present  model  using  a  non-linear  least 
squares  method.  The  characteristics  of  fabricated  and  commercial 
electrodes  are  listed  in  Table  1.  The  activation  loss  parameters, 
i/s  and  i0,  are  determined  and  listed  in  Table  2.  Both  parameters 
depend  on  Pt  loading.  The  increase  of  Pt  loading  will  cause  an 
increase  of  xj/  and  i0.  This  can  be  attributed  to  the  increase  of  the 
active  sites  for  hydrogen  adsorption. 

An  inconsistency  between  the  polarization  curves  for  the  pre¬ 
pared  electrode  with  Pt  loading  0.38  mg  cm-2  and  the  commercial 
on  E-TEK  with  Pt  loading  0.4  mg  cm-2  is  observed.  The  cell  volt¬ 
age  in  E-TEK  dropped  slightly  faster  than  the  former  electrode. 
Also,  the  value  of  i0  is  decreased  from  4.95  x  1 0-8  A  cm-2  (prepared 
electrode  with  Pt  loading  0.38  mg  cm-2)  to  2.86  x  10-8  A  cm-2  (E- 
TEK  electrode  with  Pt  loading  0.4  mg  cm-2).  This  may  be  attributed 
to  the  fabrication  process  in  the  preparation  of  electrode  with  Pt 
loading  0.38  mg  cm-2,  which  creates  a  better  particle  distribution 
of  electocatalyst.  The  localization  of  platinum  in  the  catalyst  layer 
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Table  2 

Variation  of  the  parameters  ^  and  io  with  platinum  loading  at  7=  298 1<. 


Parameters 

Fabricated  electrode  0.18  mg  Pt  cm-2 

Fabricated  electrode  0.38  mg  Ptcm-2 

E-TECK  electrode  0.4  mg  Pt  cm-2 

0.13 

0.28 

0.3 

i o  (A  cm~2) 

2.94  xlO-9 

4.95  xl0“8 

2.86  xl0“8 

can  be  improved  via  the  spraying  technique  [9].  The  values  of  i0 
estimated  are  lower  than  that  reported  by  Amphlett  et  al.  [16], 
which  are  107.6  x  10-8  A  cm-2  at  298  K.  From  these  results,  it  can 
be  concluded  that  not  only  Pt  loading  will  affect  the  value  of  i0. 

Typically  the  value  of  ^  is  in  a  very  narrow  range;  it  ranges  from 
about  0.1  to  0.5  [17].  These  minor  variations  make  experimenting 
with  different  Pt  loading  to  dramatically  change  the  voltage  simu¬ 
lated  by  the  present  model  not  a  very  productive  endeavor.  At  this 
point  it  is  useful  to  look  back  on  Eq.  (8),  as  we  can  see  from  the 
form  of  the  equation  the  only  constant  that  we  can  change  is  Iq. 
The  exchange  current  density  constant  varies  over  a  wide  range, 
and  thus  has  a  dramatic  effect  on  the  performance  of  fuel  cells  at 
low  current  densities.  Hence,  it  is  vital  to  design  fuel  cells  with  high 
exchange  current  densities. 


4.2.  Effect  of  the  temperature  on  the  performance  of  the  PEMfuel 
cell 

The  fuel  cell’s  voltage  as  a  function  of  cell’s  current  density  is 
shown  in  Fig.  3.  In  Fig.  4  cell’s  power  density  is  shown  as  a  func¬ 
tion  of  current  density  for  two  different  temperatures.  We  can  see 
from  these  figures  that  fuel  cell’s  efficiency  is  low  and  that  signif¬ 
icant  part  of  theoretical  output  voltage  is  lost  because  of  different 
losses  inside  a  cell.  We  can  also  notice  that  increase  of  fuel  cell’s 
operating  temperature  will  cause  increase  of  cell’s  output  voltage 
and  power.  If  we  increase  temperature  from  298  K  to  353  K,  fuel 
cell  voltage  would  increase  for  28%  (Fig.  3).  The  reason  for  this  is 
that  higher  temperatures  improve  mass  transfer  within  the  fuel 


Fig.  2.  Fuel  cell’s  voltage  as  a  function  of  cell’s  current  density  at  7=  298  K,  data:  (a) 
Pt  loading  =  0.1 8  mg  cm-2,  (■)  Pt  loading  =  0.38  mg  cm-2  and  (•)  E-TECK  electrode 
with  Pt  loading  =  0.4  mg  cm-2  [9];  model:  (dashed  line)  Pt  loading  =  0.38  mg  cm-2, 
(solid  line)  Pt  loading  =  0.38  mg  cm-2  and  (dotted  line)  E-TECK  electrode  with  Pt 
loading  =  0.4  mg  cm-2. 


Fig.  3.  Fuel  cell’s  voltage  as  a  function  of  cell’s  current  density  for  different  tem¬ 
peratures  at  Pt  loading  =  0.38  mg  cm-2 ;  model:  (dotted  line)  298  K,  (solid  line)  313  K 
and  (dashed  line)  353  K. 


Fig.  4.  Fuel  cell’s  power  density  as  a  function  of  cell’s  current  density  for  different 
temperatures  at  Pt  loading  =  0.38  mg  cm-2 ;  model:  (solid  line)  298  K  and  (dashed 
line)  353  K. 
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Fig.  5.  Fuel  cell’s  voltage  as  a  function  of  cell’s  current  density  for  different  fuel 
cell’s  reactants  partial  pressure  at  7=298K  and  Pt  loading  =  0.38  mg  cm-2;  model: 
(dashed  line)  Ph2  =  0.995  x  105  Pa  and  Po2  =  0.606  x  105  Pa,  (solid  line)  Ph2  = 
1.01  xlO5  PaandPo2  =  1.01  x  105  Pa. 


Fig.  6.  Fuel  cell’s  voltage  as  a  function  of  cell’s  current  density  for  different  mem¬ 
brane  humidity  at  7=298I<  and  Pt  loading  =  0.38  mg  cm-2;  model:  (solid  line) 
0  =  100%,  (dashed  line)  50%  and  (dotted  line)  0  =  30%. 

cells  and  results  in  a  net  decrease  in  cell  resistance  as  a  result  it 
improves  the  reaction  rate. 

4.3.  Effect  of  the  partial  pressures  and  membrane  humidity 

The  change  in  output  power  related  with  increase  in  partial  pres¬ 
sures  is  shown  in  Fig.  5.  We  can  see  that  power  increase  because  of 


reactant  pressure  increase  is  1 .1  %  and  is  smaller  than  change  caused 
by  temperature  variations. 

Sufficient  gas  stream  humidification  is  essential  to  PEM  fuel 
cell  operation  since  water  molecules  move  with  the  hydrogen 
ions  during  the  ion  exchange  reaction.  Insufficient  humidifica¬ 
tion  water  dehydrates  the  membrane  and  can  lead  to  cracks 
or  holes  in  the  membrane.  If  we  decrease  membrane  humid¬ 
ity  from  100%  to  30%,  fuel  cell  voltage  would  decrease  for  1.8% 
(Fig.  6).  Reducing  the  membrane  humidity  can  result  in  slightly 
slower  electrode  kinetics,  including  electrode  reaction  and  mass 
diffusion  rates.  We  can  see  from  these  simulations  that  we 
must  pay  attention  on  temperature  while  reactant  pressures  and 
membrane  humidity  are  less  significant  on  output  voltage  and 
power. 

5.  Conclusion 

Using  present  PEM  fuel  cell  model,  we  have  analyzed  the 
influence  of  fuel  cell  operating  parameters  (temperature,  partial 
pressures  and  membrane  humidity)  on  fuel  cell’s  performance.  We 
have  found  that  temperature  has  significant  influence  on  output 
voltage  and  power.  However,  the  influence  of  partial  pressures  and 
membrane  humidity  is  less  significant. 

Both  the  electron  transfer  coefficient  and  exchange  current 
density  are  platinum  loading  dependent.  The  exchange  current 
density  constant  varies  significantly  with  platinum  loading,  and 
thus  has  a  dramatic  effect  on  the  performance  of  fuel  cell  at  low 
current  densities.  This  parameter  is  vital  to  design  the  PEM  fuel 
cell. 
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